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We investigate the many-body properties of a two-dimensional electron gas constrained to the 
surface of a sphere, a system which is physically realized in for example multielectron bubbles 
in liquid helium. A second-quantization formalism, suited for the treatment of a spherical two- 
dimensional electron gas (S2DEG), is introduced. Within this formalism, the dielectric response 
properties of the S2DEG are derived, and we identify both collective excitations and a spectrum 
of single-particle excitations. We find that the single-particle excitations are constrained to a well- 
. defined region in the angular momentum - energy plane. The collective excitations differ in two 

important aspects from those of a flat 2DEG: on a sphere, the 'spherical plasmons' have a discrete 
, frequency spectrum and the lowest frequency is nonzero. 
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> ■ I. INTRODUCTION 

(N 

■ When a film of electrons on a flat helium surface reaches a critical density, the helium surface becomes unstable Q| 

, and multielectron bubbles form These multielectron bubbles are spherical cavities in the helium liquid, containing 

■ from a few to several tens of millions of electrons. The bubble radius R is determined by balancing the Coulomb 
repulsion of the electrons in the bubble with the surface tension of the helium : for N = 10000 electrons, the typical 
bubble radius is one micron, and scales as TV 2 / 3 ||. Density functional calculations indicate that the electrons 
inside the bubble are not spread out homogeneously, but instead form a thin spherical layer with thickness S <C R 
and radius as R — 5, anchored to the surface of the helium, with a binding energy of the order of sever alkelvin §|. 
The electrons are free to move in the directions tangential to the spherical helium surface, so that in effect they form 
a spherical two-dimensional electron gas (S2DEG). 

Our present analysis of the spherical two-dimensional electron gas, albeit motivated by the study of multielectron 
bubbles, is equally relevant from a fundamental point of view: it is the logical next step to take after the analysis 
of the flat two-dimensional electron gas, a topic that keeps drawing renewed attention. Furthermore, S2DEGs also 
appear in doped semiconductor particles if carriers accumulate in a surface layer in charged droplets Q, and in 
fullcrcncs, where the multipole excitation modes with angular momentum / = 1.2,3,4 of the S2DEG have recently 
been investigated ||]. 

In this paper, we introduce a second-quantization formalism suited for the description of the S2DEG (in Sec. II). 
On the basis of this formalism, the angular-momentum dependent dielectric function of the S2DEG is derived within 
the RPA framework in Sec. III. Results are given for the single-particle and collective excitations of the S2DEG and 
the static structure factor dependence on the angular momentum quantum number in Sec. IV. 
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II. HAMILTONIAN OF THE S2DEG 



Decompositions in plane waves <fik(r) oc e lk r are well suited to handle the many-body problem in flat space, but are 
not appropriate to describe the many-body problem on a spherical surface. The non-interacting single-particle wave 
functions suitable for the description of the S2DEG are spherical harmonics, ipi tTn (fl) = YJ. m (f2), where f2 represents 
the spherical angles Q — {0, </>}, and l,m represent the angular momentum quantum numbers. We introduce the 
creation operator cf m which creates an electron in the state v?;, m (Q) on the sphere if no electron is present in that 
state, and the annihilation operator Q. m which destroys an electron in the single-particle state tpi tm {VL). 



1 



A. Kinetic energy 



Within second quantization theory (see e.g. pj), the kinetic energy operator can be written as a function of . 
and ci m as follows 
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(1) 



In this expression, to c is the electron mass, R is the radius of the spherical electron gas, Aq is the angle-dependent 
part of the Lagrangian in spherical coordinates, and f dfl = ^ d(j> d9 sin 9. Expression (|]) shows that the energy 
of a free electron in an angular momentum state {I, to} is 



l,m 



h 2 l(l + 1) 

2m r R 2 



(2) 



If more than one electron is present, they will occupy the lowest-energy states up to the Fermi energy Ep characterized 
by a Fermi angular momentum Lp, so that N — 2(Lp+ l) 2 and Ep = H 2 Lp(Lp + 1) / (2m c R 2 ) . Since in multiclcctron 
bubbles, the bubble radius R cx N 2 ^ 3 p|, the Fermi energy scales as Ep cx N~ li/3 so that adding electrons to the 
multielectron bubbles lowers the Fermi energy. 



B. Density 



The surface density on the sphere, defined as the expectation value of the number of electrons in an infinitesimal 
spherical angle around il, is 



oo I oo V 

h e (n) = j2 E E E Y l : m (m',rrAn)c+ m c l , 
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The decomposition in spherical harmonics of this density operator is given by 
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«e(^) = E E pS^ylm, 
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p e (l,m) — / dtl h e (rt)Yi, m (tt). 
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Using the parity property of spherical harmonics, Yi :m (Jl) = (— l) m Y*_ m (n), it is easily shown that p e (l,m) 
(— l) m /5 c (7, — m). Substituting (0) in (Eh, we find for the spherical components of the density operator: 
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The addition rules for spherical harmonics allow to perform the integral of the product of spherical harmonics in 

' ,+ '' 1 1 {21 + 1){2V + 1) 
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In this expression, (I, m; I', m'\L, M) is the Clebsch-Gordan coefficient for combining the angular momenta {/, m} and 
{7', to'} into {L,M}. To simplify this expression, we introduce the following notation: 
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i+V 

E 

L=\l-l'\ M=-L 



^,m)m',m')= E E \/ ~ + { 2L+'i) ^ {l > 0; l '' 0|£ ' 0) {l ' m; ^ W/|L ' M) ^M- (8) 



The operator m \g,/j, m n creates an electron in a state which results from the combination of two angular momentum 
states characterized by the quantum numbers {I, m} and {l',m'}, respectively. With this notation, the spherical 
components of the density operator can be written as 

oo I 

p e (l,m) = Y^ E £ t™)m',™>fi',m', (9) 

!'=0m'=-l' 

which highlights the analogy with the definition of the Fourier transformed density operator of the flat two-dimensional 
electron gas (2DEG). 



C. Coulomb interaction on the sphere 



The Coulomb potential energy term for N electrons, in the Hamiltonian is, 

N N 



-^coulomb — ~ ^ ^ "j TZ ^-7 j (10) 



where f j represents the position operator of electron j and £o is the vacuum permittivity. The Coulomb potential on 
a sphere with radius R can be straightforwardly expanded in spherical harmonics: 

1 N N 2 00 1 a 

H coulomh = -Y: E i^E E uhYt^Y^cij,), (ii) 

J=lj'^j'=l i=0 m=-/ 

with fij the (spherical angle) position operator of electron j on the sphere. In second quantization this becomes 

^coulomb = ^ ^ XX - !)" 1 U (0 ^(Z 1 ,mi)(8l(i,-m)^(Z2,m2)®(i,m)^ 2 ' m2 ^ 1 ' mi ' ( 12 ^ 
Ji.mi i 2 ,rri2 '' m 

With E(,m = Ej=0 E!n=-( and 

" 5SSTT < 13 » 

This term in the Hamiltonian describes the interaction process, shown in Fig. 1. The initial state consists of two 
electrons in angular momentum states {/i, m{\ and {I2, WI2}. In the final state one electron is in an angular momentum 
state resulting from adding {I, m} to its original angular momentum, and the other electron is in a state resulting 
from subtracting {l,m} from its original angular momentum. Hence, the Coulomb interaction on a spherical surface 
can be described as the exchange of a virtual photon with a given angular momentum {l,m}, with an amplitude for 
this process of v(l). The total Hamiltonian for the S2DEG is 

1 — 771— — I 

+ E E £(-ir«(o 

(14) 

h,mi l 2 ,m 2 L ' m 
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III. RESPONSE PROPERTIES OF THE S2DEG 



When an external field V QX t(Q.;uj) = V ex t(l, m; oo)Yi im (il) that couples to the electron density is applied, an 

induced density p- ln< i(^l,uj) = J2i m P'md(l, m\ uj)Yi im (n) is generated. Within linear response theory, the spherical 
components of this induced density can be written as 



pind(l,m;cj) = ^ e7A {l,m-,u)T> R (l,m]oS), 



(15) 



where T> r(1 , m; lj) is the retarded density-density Green's function fl2]| . If the external potential V cxt (l , m; u>) is a 
Coulomb potential created by an external charge p ex t, a dielectric function depending on the angular momentum can 
be introduced to describe the screening of this external Coulomb potential: 



e(l, m; u>) 



p ind (l,m;uj) 
p cxt {l,m;uj) 
1 . v(l) 



= 1 



T>fi(l, m; u>). 
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The retarded Green's function T>r(1,iti;lu) is derived from the density-density Green's function 



T>(l,m;t) 
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(16) 
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where the expectation value is taken with respect to the many-body ground state |Wq), T is the time-ordering operator, 
and 



Cl,m(t) 



_ iHt/H 



-iHt/H 



(19) 



Within the second quantization formalism for the S2DEG, T> r(1 , m; lj) is evaluated with standard Green's function 
techniques. In the resulting Feynman graphs (the 'polarization bubble' graphs, |]l]|), the wave numbers are replaced 
by angular momenta as illustrated in Fig 1. To lowest order in the interaction amplitude v(l), we find 
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with 77 an infinitesimal number, and Ei jlri the energy of a free electron in angular momentum state {/, m} on the 
sphere and 



n(l, m) 



(21) 



the occupation number of state {l,m}. This result, derived with the formalism described in Sec. II, is in agreement 
with the result of Inaoka (^] for the susceptibility of the S2DEG, which was used to investigate the I — 1,2,3,4 
multipole modes of the S2DEG |l. In Sec. IV we extend this investigation to include excitations with large angular 
momentum I (so that l/Lp > 1), keeping the regime studied by Inaoka as a limiting case. By investigating a much 
broader portion of the excitation spectrum, we identify the nature of the excitations and expose novel properties of 
those excitations. 

The contribution of lowest-order in v(l) to the retarded density-density Green's function determines the Hartree- 
Fock approximation (HF) to the dielectric function (|l7j), given here for reference: 



£ H fG,»i; oj) 



1 1 "(0^(0)/, \ 

1 + — Vy(l,m;oj) 



(22) 



It also permits the calculation of the random-phase approximation (RPA) to the dielectric function, through a Dyson 
series for 2?^ PA : 
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T)(0) 

vf A = : r ^„ (or (23) 



i - K0/?TO 

so that 

eRPA (/,m;c) = l-^< ) (/,m; W ). (24) 

Note that the spin degree of freedom is not explicitly taken into account in the expressions above. The interparticle 
potential and the kinetic energy are considered to be spin-independent in the present treatment. The central quantity, 
from which properties of the S2DEG are derived in the current treatment, is the 'polarization bubble' diagram leading 
to T>jp . For spin-independent interaction potentials, the spin degree of freedom only leads to a degeneracy prefactor 

-,(0) 



T> R : the degeneracy of the unperturbed levels is doubled from 21 + 1 to 2(21 + 1). 



IV. RESULTS AND DISCUSSION 

A. Single-particle excitations 

In this section we will use units such that m c = H = R = 1. The dielectric function can be clarified by considering 
Plemelj's rule. The imaginary part of the RPA dielectric function becomes 

/+ '' (21 + l)(2l' + 1) 

Im e-rp A {1, m; uj) = nv (I) ^2 ^2 — An(2L + 1) — n ( l ' m ')i 1 - n ( L > m + m ')] ( 25 ) 

I'm' L=\l-l'\ ^ ' 

x 0; I', 0\L, 0)\ 2 \(l, m; I', m'\L, m + m')\ 2 5[w - {E L<m+m , - E v , m ,)]. 
The dynamic structure factor of the S2DEG is related to the dielectric function by 



S(l,m;u) = TTrlm 

irv{l) 



e(l, m; w) 
1 lm[e(l,m;u))] 
™(0 {Re[e(/, m;w)]} 2 + {Im[e(i, m;u)]} 



(26) 
(27) 



and can be interpreted as the probability that an excitation with given angular momentum quantum number {l,m} 
and energy hu can be created. S(l, m;uj) is a quantity accessible to experiment, in particular scattering experiments. 
From expression ( |25| ) we conclude that ImeRPA(/,m;o;) is zero unless adding the angular momentum {l,m} and the 
energy Hui to the ground state of the S2DEG can excite a single electron from an occupied state [n(l', m') = 1] into an 
unoccupied state [n(L,m + m') — 0]. These are the single-particle excitations (called the "Landau continuum" in the 
case of the flat 2DEG). Fig. 2 illustrates this concept in relation to the S2DEG. Three regimes can be distinguished: 

• Case 1: I < Lp. In the left panel of Fig. 2, we show the case of angular momentum {I, m} = {2, 1} imparted 
on the system. The filled disks represent the occupied states, the open squares are allowed final states which 
can be reached by adding the angular momentum {I, m} to the angular momentum of an electron in the Fermi 
sphere, and the open circles are forbidden final states as they are occupied and excluded by Fermi statistics. In 
this example the Fermi sea is filled up to Lp = 3. Not all final states are accessible, since to excite an electron 
it has to go into a final state which is unoccupied (the open squares) . This means in practice that the electrons 
which can participate in creating an excitation of angular momentum {I, m} are the ones close to the Fermi 
level (in fact, those from level Lp — I up to hp). 

• Case 2: 2Lp > I > Lp. This is shown in the middle panel of Fig. 2 depicting single-particle excitations which 
have angular momentum {4,-3}. Now all the electrons in the Fermi sea can participate, but not all resulting 
final states are unoccupied. 

• Case 3: I > 2Lp. An example is shown in the right graph of Fig. 2, for single-particle excitations of angular 
momentum {8,1}. Now not only can all electrons participate in the process, but also all possible final states 
are unoccupied and thus accessible. 
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The maximum energy difference between initiai and final single-particle states in Fig. 2 is 



2m R 2 

The smallest energy difference in case 3 (I > 2L F ) is 



bW(0 = + l )i L F + 1 + 1)- (L F )(L F + 1)] 

"~ [l 2 + 1{2L F + 1)] . (28) 



Wmi»(0 = 2^jl?2 [0 - £ f )(' - ^ + 1) - (if)(if + 1)] 

[l 2 - l(2L F - 1) - 2i F ] . (29) 



2m c R 2 



The two frequencies u) m i n (l) and w max (Z) demarcate a region in the frequency versus angular momentum plane in 
which the single-particle excitations lie. Fig. 3 shows the location of the excitations of the spherical two-dimensional 
electron gas in the angular momentum I vs. frequency u) plane. The {I, to} values corresponding to a single-particle 
excitation for an S2DEG with L F = 10 are shown in Fig. 3 as a markers within the limiting frequencies w m j n and 
w max of the Landau continuum region shown by full curves. 

B. Plasmons in a spherical 2D electron gas 

In addition to the single-particle excitations in the S2DEG, collective excitations are also possible. "Collective 
excitations" are defined as poles of the density-density Green's function T>, whereas single-particle excitations are 



defined as poles of the single-particle Green's function |12|. This means that collective excitations appear when 



e{l,m\bj) = 0. In the flat 2DEG, these collective excitations are called plasmon modes. We will use the same 
terminology for the S2DEG (occasionally using 'spherical plasmons' when a distinction is needed). 

For the collective mode, the imparted angular momentum is shared between all the particles. This means that in a 
collective mode with angular momentum {I, m} the entire spherical shell of electrons will oscillate with an amplitude 
proportional to Yj iTO (il). The frequency of the collective modes (the plasma frequency uj p \) will depend on the angular 
momentum of the mode m}. The plasma frequency u> p i(l, m) of the S2DEG can be found by solving 
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-Re[T> Ri0 (l,m;u p i)]=Q _ ^q) 



1 2ehR 21 
lm.[T>R ! o(l,m;u}p{)] = 

For the flat electron gas the last condition means that the plasma branch lies outside the region of single-particle 
excitations discussed in the previous subsection. In Fig. 3, the spherical plasmon excitations are shown as full 
diamonds. Fig. 4 illustrates how the plasmon branch depends on the number of electrons in the S2DEG. The radius 
of the S2DEG has been chosen equal to the equilibrium radius of a multielectron bubble with a given number of 
electrons pi] . The bubble radius is found by balancing the helium surface tension and the Coulomb repulsion as 
discussed in the introduction. Fig 5. illustrates how the plasmon branch depends on the radius of the S2DEG, for 
a fixed number of electrons. We find that as the radius is decreased (i.e. the multielectron bubble is compressed), 
the plasmon branch lies relatively closer to the upper frequency of the single-particle excitation region (as shown 
in the inset of Fig. 5). However, the energy scale, set by H 2 /(m e R 2 ), increases with decreasing radius so that the 
absolute value of the plasmon frequencies increases with decreasing radius. The plasmon frequency is independent of 
the projection of the angular momentum m. The plasmon branch of the S2DEG strongly differs from the plasmon 
branch of 2DEGs in two important aspects: it is discrete and the lowest accessible plasmon frequency is larger than 
zero despite the acoustic nature of the plasmon branch. 

C. Sum rule and spectral weights 

Sum rules prove to be a remarkably useful tool in the analysis of spectra, both experimentally and theoretically. 
The zeroth moment sum rule for the dynamic structure factor defines the static structure factor: S(l,m) : 



G 



S(l, to) 



S(l, to; u>)du> 
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j Im 


e(Z, to; uS) 





dio. 



(31) 



Since the dynamic structure factor only depends on the magnitude of the angular momentum I and not on to, this will 
also be true for the static structure factor. Fig. 6 shows the static structure factor as a function of angular momentum 
I. The open squares are the result in the Hartree-Fock approximation [using £hf, expression (p2|)], and the full circles 
are the result in the RPA approximation [using £rpa, expression (p4|)]. As in the case of the flat 2DEG, the HF 
structure factor S*hf(0 is linear for small l/Lp, whereas the RPA structure factor approaches zero more rapidly with 
decreasing I. In the inset of Fig. 6, the first frequency moment of the dynamic structure factor, 



MO) 



(32) 



is shown, again for both HF and RPA approximations. Inaoka || derived a sum rule for the first frequency moment 
of the angular momentum dependent dynamic structure factor. In units h — m e = R, this is 



m = 



(33) 



From Fig. 6 it is clear that the HF dynamic structure factor obeys this sum rule. The RPA dynamic structure factor 
can be written as a sum of a contribution from the plasmon mode and the contribution Spoilt from single-particle 
excitations: 



S rpa (1;cj) = A(1)S[lj - w p i(i)] + S cont (l;u), (34) 

with A(l) the spectral weight of the plasmon branch. The inset of Fig 6. shows J uiS con t(l,uj)dui. From the deficit of 
the RPA dynamic structure factor without the plasmon mode, shown in the inset of Fig. 6, the spectral weight of the 
plasmon mode can be derived: 

oo 

A(l)w pl (l) = l !±}± _ f uS cont (l,u;)dLJ. (35) 
o 

The spectral weight of the plasmon mode, along with the RPA dynamic structure factor in the region of single-particle 
excitations, is shown in Fig. 7 for Lp — 20. The RPA dynamic structure factor of the S2DEG is shown color-coded 
and as a function of the angular momentum I and the energy hw. In white regions the RPA dynamic structure 
factor is zero. In colored regions, the RPA dynamic structure factor differs from zero. This figure is the complement 
of Fig. 4: whereas Fig. 4 shows the location and type of the possible excitations of the S2DEG as a function of 
angular momentum and energy, this figure shows their corresponding spectral weight as expressed by the dynamic 
structure factor. At small angular momentum (l/Lp < 0.5), the plasmon branch carries the most spectral weight, 
and the region of single-particle excitations only has a small fraction of the total spectral weight. Within the region 
of single-particle excitations, there is a local maximum in the dynamic structure factor around I = Lp- 



D. Discussion 



The dynamic structure factor S(l, to; u>) is the spectral function of the density-density Green's function and as such 
can be interpreted as a probability function (see e.g. Jr|, p. 153). S(l, m; u)du> is the probability that an excitation of 
the S2DEG has angular momentum quantum numbers {I, to} and energy between hu and h(u> + du)). This is analogous 
to the interpretation of the more familiar dynamic structure factor [S*(k, lu)] expressed as a function of wave number, 
which is the probability that an excitation has momentum Kk and energy hu. 

A 'natural' experiment to determine the dynamic structure factor directly is an inelastic scattering experiment. In 
particular the differential cross section for a probe particle of mass M (such as an incident ion) and initial momentum 
Pi to be scattered to a final state pf with energy transfer between huj and K(u) + duj) to the scattering system is [ ^2[ : 

^ =^WS(*,«). (36) 



dOduj 8tt 3 pi 
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In this expression, 9 is the scattering angle between pt and p;, \Vk\ is the Fourier transform of the interaction 
potential between the probe particle, and the energy transfer to is related to the momentum transfer k by fko — 
ft 2 k.pi/m + (hk) 2 /(2m). The relation between the Fourier decomposition and the spherical decomposition of the 
density is 

N N 

p e (k) = (e lk ' ?J ) =X)EVM2/ + iyjl(*i2) (Yi,o@j)) 

3=1 3=1 I 

= VMM + l)i l 3i(kR)p e (l, 0) (37) 
i 

where the z-axis has been taken along the direction of the wave vector k, and ji(x) is the spherical Bessel function 
of the first kind, and 9j is the angle between k and the position operator of electron j, Tj. The dynamic structure 
factors are related by 

S(k, u) = | VMM + ( k R) | 2 S(l, 0; u) (38) 
l 

which can be derived by writing the dynamic structure factor as a density-density correlation function in the Lehmann 
representation and substituting expression ( |37| ) for p e (k). 

Note that the dependence of the angular wave functions on the polar angle 6 becomes semiclassical when the 
quantum number I is large. In particular, the WKB expression for the angular function becomes Jl4[ 

Yu t *M + w + *f4 (39) 

7T V sin 9 

This semiclassical approximation will be valid in the region where 91 3> 1 and {it — 9)1 3> 1, which holds at large I 
for almost all values of 9 except those close to the 'poles' 9 = and 9 = tt. The semiclassical approximation ( |39| ) 
corresponds to the limit of a plane wave, where the wave length is given by R(2n/l). As such, this approximation 
leads primarily to the results of a flat 2D electron gas where hl/R plays the role of the electron momentum. This is 
valid in the limit that R is large and the ratio l/R remains finite, although even for small R, this approximation will 
still hold at large I. However, it does not reproduce some of the interesting low angular momentum results such as a 
lower boundary to the plasmon frequencies. 



V. CONCLUSIONS 



In this paper, we have studied the properties of the two-dimensional electron gas in a distinct geometry which is 
recently gaining interest in connection to multielectron bubbles, namely electrons confined to a spherical surface. For 
this purpose, we set up a second-quantization formulation based on the spherical harmonics as single-electron building 
blocks for the many-electron theory, introduced in Sec. II. Within this formalism, the dynamic structure factor is 
derived in the RPA framework and both the single-particle excitations and collective excitations are analyzed for 
S2DEGs with Lf up to 40 and for angular momenta I up to 3Lf ■ 

The dynamic structure factor in the RPA approximation reveals collective excitations ('spherical plasmon modes'), 
which differ from the collective modes of the flat 2DEG in that the spherical plasmon modes of the S2DEG are 
discrete in frequency, and the smallest spherical plasmon frequency is larger than zero. For l/Lp < 0.5, we find 
that the spherical plasmon mode carries the main fraction of the spectral weight. The single-particle excitations 
are confined to a region determined by u> m m(l) and w max (Z) given by expressions ( p9| ) and (p8|). In the semiclassical 
approximation, the role of the electron momentum hk in the flat 2DEG is played by hi/ R in the spherical electron gas, 
and the energy scale is set by h 2 j (m e R 2 ) , which opens interesting prospects for multielectron bubbles since the radius 
can be varied by cither changing the number of electrons or -independently- the pressure. For typical multielectron 
bubbles with TV ^ 10 4 , the plasmon modes of the S2DEG lie in the far-infrared, and their frequency increases with 
decreasing number of electrons in the multielectron bubble or with decreasing radius of the bubble, such that these 
novel collective modes may be detectable in forthcoming experiments on stabilized multielectron bubbles p5[ . 
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FIGURE CAPTIONS 



FIG. 1. This figure shows two Feynman diagrams for equivalent processes, namely the exchange of a virtual photon 
between two electrons. The panel at the left shows the diagram for the case of a flat 2DEG and the one to the right 
is for the case of a S2DEG. In the flat 2DEG, wave numbers are good quantum numbers, in the S2DEG they have to 
be replaced by angular momentum quantum numbers. Vector addition of momenta has to be replaced by addition of 
angular momenta. 

FIG 2. This figure shows diagrams representing the angular momentum states of an electron on a spherical surface. 
In the three panels, electrons are occupying the single-particle states up to a Fermi-level angular momentum Lp = 3. 
This Fermi sea of occupied states is shown as filled circles. The diagrams also show the final states which can occur 
from combining the angular momentum of any of the electrons in the Fermi sea with the angular momentum of an 
excitation (with {l,m} — {2, 1}, {4, — 3}, {8, 1} for the left, middle and right panels, respectively). The final states 
which are unoccupied, and thus allowed, are indicated by hollow squares. The final states which cannot be achieved 
because there is already an electron present are indicated by hollow circles. 

FIG. 3. In the many-body S2DEG, two types of excitations can be distinguished: single particle excitations (which 
appear for Im[e] ^ 0) and collective excitations (which appear when Re[e] = Im[e] = 0). The excitations of the 
S2DEG can be characterized by their angular momentum I and energy Hlo. This figure shows the location of the 
excitations in the frequency (energy) vs. angular momentum plane. The single particle excitations are discrete, and 
all lie in a region demarcated by the black curves, representing u) m i n (l) and w max (0 given by expressions (p9|) , (p8|) . 
The collective excitations are shown as filled diamonds, lying above the w max (Z) curve. 

FIG. 4. The collective excitations of the S2DEG are the analogue of plasmons in the flat 2DEG and appear at 
frequencies u> and angular momenta I such that Re[e(/,tj)] = Im[e(Z,w)] = 0. For S2DEGs with Lp = 10,20,40, 
these collective excitations are shown as filled squares, circles and triangles, respectively. The dashed curves denote 
the upper frequency of the region of single-particle excitations (given by expression (|28|)). As more electrons are put 
into the multielectron bubble, the frequency of the collective excitations decrease. As the Fermi level is increased, 
Lp — > oo, the spherical plasmon excitations form a continuous curve and the lowest spherical plasmon frequency 
u) p \(l = 1) approaches zero. 

FIG. 5. In a multielectron bubble (MEB), the radius of the spherical electron gas can be varied independently of 
the number of electrons by applying pressure. This figure, which complements Fig. 4, shows the effect of changing 
the MEB radius on the frequencies of the collective excitations (the 'spherical plasmons') as a function of the angular 
momentum. The inset shows the same results, with the frequency rescaled to the natural bubble frequency h/(m c R 2 ). 

FIG. 6. The static structure factor S(l) for a spherical two-dimensional electron gas (S2DEG) with Lp = 20 is 
shown as a function of I, both in the RPA and in the Hartree-Fock approximations. Note that in the Hartree-Fock 
approximation, S(l) is linear in I at small I, whereas in the RPA approximation it goes to zero faster than linearly. 
In the inset, the sum rule ([33]) for the first frequency moment of the dynamic structure factor of the S2DEG || is 
checked. The Hartree-Fock result complies perfectly to the sum rule. The RPA result without the plasmon branch 
does not have enough spectral weight to satisfy the sum rule for l/Lp < 2, which indicates that in this region a 
substantial part of the spectral weight lies with the plasmon mode - the strength of the plasmon branch is derived 
from this result. 

FIG. 7. The RPA dynamic structure factor S(1,uj) for a spherical two-dimensional electron gas (S2DEG) with 
Lp = 20 is shown here in units h = m c = R = 1 as a function of the angular momentum quantum number I and their 
energy u. At zero temperature, the RPA dynamic structure factor of the S2DEG consists of a region of single-particle 
excitations and a set of collective modes (colored disks). This figure complements Fig. 3, which shows only the 
locations of the excitations, whereas this figure shows their spectral weights. 
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